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1.0  INTRODUCTION 

In  high  enthalpy  gaseous  flows  associating  high  velocities  and/or  high  temperatures,  physical  and 


chemical  processes  such  as  vibrational  excitation,  dissociation,  ionisation  and  various  reactions,  can  take 
place.  The  characteristic  times  of  these  processes  have  often  the  same  order  of  magnitude  as  the 
“mechanical”  or  aerodynamic  characteristic  times,  so  that  these  flows  constitute  typical  non-equilibrium 
media. 

The  best  way  for  analysing  these  reactive  flows  in  continuous  or  “collisional”  regime  consists  in  using  a 
statistical  approach  by  considering  the  macroscopic  quantities  as  local  averages  of  various  properties  of 
elementary  particles  (molecules,  atoms,  ions,...)  and  by  taking  into  account  their  interactions  resulting 
from  their  “collisions”.  Thus,  the  Boltzmann  equation  seems  to  be  an  appropriate  tool  for  the  description 
of  these  flows. 

Among  the  methods  used  for  solving  the  Boltzmann  equation,  the  Chapman-Enskog  method,  consisting  in 
expanding  the  distribution  function  in  a  series  of  a  “small  parameter”  equal  to  the  ratio  of  the 

characteristic  time  between  collisions  to  a  reference  flow  time  has  known  a  great  success.  However,  in  the 
past,  it  has  been  generally  limited  to  the  case  where  only  one  single  type  of  collision  is  present  in  the 
medium.  Thus,  at  relatively  low  temperature,  when  the  elastic  collisions  are  “dominant”,  the  behaviour  of 
the  system  is  correctly  described  by  Navier-Stokes  equations  in  which  transport  terms  are  satisfyingly 
calculated  by  the  Chapman-Enskog  method  [1],  [2],  [3]. 

Now,  when  physical  and  chemical  processes  take  place,  the  ratio  of  their  characteristic  times  sn  (inelastic 
and/or  reactive  collisions)  to  the  reference  flow  time  can  take  any  value.  The  problem  first  is  to  compare 
Sj  and  su  and  to  insert  the  terms  of  physical  and  chemical  production  in  the  hierarchy  imposed  by  the 

expansion  in  a  series  of  £j  and  then  to  compute  the  modifications  brought  to  the  Navier-Stokes  system 
and  to  the  transport  terms  [4],  [5]. 

It  is  obvious  that  a  great  variety  of  non-equilibrium  situations  may  exist  due  to  the  numerous  possible 
multi-scale  physical  and  chemical  processes.  However,  it  is  imperative  that  the  number  of  collision  types 
should  be  restricted  to  only  two  (collisions  I  and  II),  in  order  to  avoid  an  expansion  beyond  the  first  two 
terms  in  the  Chapman-Enskog  expansion  (Navier-Stokes  level)  [6]. 

The  strict  application  of  the  Chapman-Enskog  method  to  non-equilibrium  situations  is  first  presented:  it 
leads  to  two  main  general  approaches  called  WNE  and  SNE  methods  depending  on  the  degree  of  non¬ 
equilibrium  considered  :  WNE  for  “weak”  non-equilibrium  situations  and  SNE  for  “strong”  ones  [7].  In 
the  application  of  these  methods  to  concrete  cases,  it  is  obvious  that  they  do  not  match  in  the  intermediate 
situations.  That  is  why,  a  generalized  Chapman-Enskog  method  (GCE),  capable  of  realizing  this  matching 
is  finally  presented  and  developed  [8],  [9]. 

On  the  other  hand,  as  it  is  impossible,  in  the  framework  of  this  lecture  to  analyse  all  possible  situations, 
two  typical  situations  are  examined  hereafter:  first  the  case  of  vibrational  non-equilibrium  in  pure  gases 
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and  gas  mixtures  and  second  the  case  of  dissociating  pure  gases.  A  few  results  are  also  given  for  more 
complex  systems  [10],  [11]. 

At  last,  applications  of  the  results  obtained  from  the  GCE  method  are  presented  for  hypersonic  flows  in 
which  vibrational  and  chemical  non-equilibriums  simultaneously  take  place:  they  concern  the  kinetic  rate 
constants,  the  transport  coefficients  and  more  generally  the  characteristics  of  hypersonic  flows  such  as 
shock  waves,  boundary  layers  and  flows  around  bodies  [12],  [13],  [14]. 


2.0  GENERALITIES  ON  NON-EQUILIBRIUM  FLOW  REGIMES 


2.1  The  Boltzmann  Equation  in  Reactive  Systems 

In  a  flow  of  gas  mixture,  the  Boltzmann  equation  governing  the  distribution  function  f.  of  the 
particles  of  the  component  p  in  the  quantum  state  i  can  be  simply  written 

dfip 


dt 


=  J 


(1) 


where  the  quantum  level  i  includes  a  rotational  level  ir ,  a  vibrational  level  iv  and,  eventually  an 
electronic  level  ie :  i  =  (ir,  iv ,  ie ) 

The  collisional  term  J  includes  all  types  of  collision,  elastic,  inelastic  and  reactive,  that  is 

J  —  J  j  +  J  R  +  Jy  +  J Q  H“  ...  ,  (2) 

where  JT  represents  the  collisional  balance  of  particles  undergoing  a  velocity  change  only  (elastic 
collisions),  JR  and  Jv  the  balance  of  collisions  with  rotational  and  vibrational  change  respectively 
(inelastic  collisions),  J c  the  balance  of  collisions  with  species  change  (reactive  collisions),  etc. . . 


These  collisions  may  be  subdivided  in  more  specific  types,  such  as  VV  collisions  (both  interacting 
molecules  change  their  vibrational  level)  ,  TV  (only  one  molecule  change  its  level)  or  resonant  collisions 
(they  interchange  their  level),  etc. . .. 

These  collisions  of  course  have  different  probabilities  to  occur,  so  that  they  have  different  characteristic 
times  (probable  times  between  specific  collisions),  rr,  tr ,  zv ,  rw,  rc ,...  which  may  differ  by  an  order 
of  magnitude.  Thus,  two  examples  are  given  in  Figs.l  and  2  for  02  and  N2  respectively. 

In  these  conditions,  we  may  write  the  Boltzmann  equation  (1)  in  the  following  dimensionless  form: 


dfr 


V 


1 


1 


1 


1 


J j  +  J ^  +  Jy  +  J q 

dt  Sj  £ p  Sy  £ q 


(3) 


or  any  other  form  adapted  to  the  specific  case  considered. 

/  £  =  Tr/ 

'9  ’  R  /O  ’ 


In  Eq.(3),  we  have:  sT  =  Tt/  ,  =  Tr/  ,  £T/  =  Ty/ 


where  6  represents  the  “aerodynamic”  reference  time,  depending  on  the  problem  under  study  and  that 
may  take  considerably  different  values. 


3-2 


RTO-EN-AVT-1 62 


Kinetic  Theory  of  Reactive  Molecular  Gases 


As  stated  in  the  introduction,  we  will  consider  only  two  collisional  characteristic  time  scales  and  £n  , 

so  that  all  collisions  are  classified  in  two  groups,  for  example  elastic  collisions  in  group  I  and  inelastic  in 
group  II;  or  TRV  collisions  in  group  I  and  reactive  in  group  II,  etc...  Thus,  we  will  simply  write  Eq.3  in 
the  following  form: 


dt* 


-L  *  J-  r* 

- - J1  + - JL 


(4) 


rp  ( atm.s ) 


Figure  1:  Translation,  rotation,  vibration,  dissociation  characteristic  times  (Nitrogen) 


6  8  10  12  14  T  (10  3AT) 

Figure  2:  Translation,  rotation,  vibration,  dissociation  characteristic  times  (Oxygen) 
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2.2  Chapman-Enskog  expansion:  Flow  regimes 

Considering  flows  in  “collisional  regime”,  we  may  assume  that  the  collisions  I  are  dominant  (the 
most  probable),  so  that  Sj  □  1,  (r,  □  6)  An  this  case,  we  can  classically  expand  fip  in  a  series  of 

*,=*.  that  is  &=  ft +*%+'■'  , 

or,  stopping  the  expansion  at  the  first  order: 

/; =4°'  (!+««;) .  <5> 

where  fip  is  the  zero  order  distribution  function  and  f'p  the  first  order  one,  with  fip  =  f°<pip  and 

1. 

Now,  we  have  to  substitute  the  expansion  (5)  in  Eq.(4)  and  to  introduce  the  collisional  term  II 


into  the  expansion.  In  order  to  do  that,  there  are  three  possibilities:  l)enU  1,  that  is eu  □ 
2)  snU  1 ,  that  is  sn  L  e ,  or  3)  e„Ul. 

In  the  first  case,  we  obtain  the  following  system  giving  successively  f?  and  (pip : 

J°  =  0 

C° 


i 


dt 


=  J\ 


(6) 


Thus,  the  collisions  II  have  no  influence  on  the  zero  and  first  order  distribution  function.  This  corresponds 
to  the  “frozen  case”  for  the  collisions  II. 


On  the  contrary,  in  the  case  2),  the  collisions  II  play  the  same  role  as  the  collisions  I  in  the  determination 
of  the  distribution  function,  since  we  have  the  following  system  to  solve: 


(J,+J«f=  0 

<=(J  +J  y 
dt  ' ' 

In  the  intermediate  case  3),  we  have  the  following  system: 

J°j  =  0 

dfl 


dt 


.  =  Jl  +  J° 

J  I  T  J  u 


(7) 


(8) 


Thus,  we  can  dismiss  the  frozen  case  which  brings  nothing  new  in  comparison  with  the  classical  CE 
method  used  when  only  one  type  of  collision  is  considered. 


Considering  now  the  case  2),  we  easily  find  that,  at  the  zero  order  of  the  system  (7),  the  Euler  equations 
govern  the  system  with  an  equilibrium  distribution  (levels  and  species)  for  the  zero  order  distribution 

function  .  On  the  contrary,  in  the  case  3),  at  zero  order  of  the  system  (8),  the  Euler  equations  must  be 

closed  by  kinetic  equations  (species)  and/or  relaxation  equations  giving  the  evolution  of  “populations”  due 
to  collisions  II  in  the  time  scale  6 :  this  corresponds  to  a  non-equilibrium  situation. 
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At  first  order,  the  linear  integral -differential  equations  allow  to  obtain  the  perturbation  (pip  and  the  Navier- 

Stokes  equations  govern  the  system,  but  kinetic  and/or  relaxation  equations  are  still  necessary  in  the  case 
3)  and,  moreover,  the  collisions  I  and  II  appear  in  a  different  way  in  cases  2)  and  3). 

It  is  also  to  be  noted  that,  whatever  the  form  of  fip  is,  in  equilibrium  or  not,  the  first  order  solution  always 

exhibits  a  non-equilibrium  which  remains  small  because  of  the  linearisation.  Therefore,  the  solution 
obtained  from  Eqs.(7)  corresponds  to  a  weak  non-equilibrium  case  called  WNE  regime  and  the  solution 
corresponding  to  Eqs.(8)  is  denoted  SNE  regime  (strong  non-equilibrium). 

The  cases  presented  above  may  hide  complex  situations:  first,  the  indices  I  and  II  represent  various  types 
of  collisions,  either  simple  (T,  R,  V,  C),  complex  (TR,  TRV,...),  or  particularized  (VV,  resonant...)  ;  thus 
the  specific  situations  are  numerous. 

Furthermore,  each  type  of  regime  may  apply  to  one  or  several  types  of  collisions.  For  example,  we  call 
(WNE)I+(SNE)II  a  weak  non-equilibrium  regime  for  the  collisions  I  and  a  strong  non-equilibrium  for  the 
collisions  II.  If,  as  a  standard  example,  we  consider  the  case  of  a  pure  gas  in  a  strong  vibrational  non¬ 
equilibrium,  we  call  this  case  (WNE)TR+(SNE)V.  In  the  same  way,  if  the  chemical  non-equilibrium  is 
significant  at  the  zero  order  and  if  the  vibrational  non-equilibrium  is  “less”  important,  we  call  this  regime 
(WNE)TRv+(SNE)c,  or  simply  (WNE)v+(SNE)c,  etc. .  .Simple  typical  examples  are  presented  below. 


3.0  VIBRATIONALLY  RELAXING  GASES 
3.1  Zero  Order  Solutions 

Pure  diatomic  gas  flows  are  first  considered.  The  equilibrium  zero  order  solution  (Maxwell- 
Boltzmann)  is  well  known  and  is  obtained  from  the  equation:  JRRV  =  0 .  Thus,  the  “normal” 
solution  is 


r  0 

fi  =n 


f  m  ^ 

exp 

r  mu2  ^ 

8ir  exPl 

{-sJkT)  exp( 

v2  7tkT  j 

{  2kT ) 

<2* 

Qv 


(9) 


with  usual  notations:  n  number  density,  m  molecular  mass,  T  temperature,  sj  ,  st  rotational  and 

vibrational  energies  per  molecule,  QR(T),  QV(T) rotational  and  vibrational  partition  functions,  g 

statistical  weight.  One  single  temperature  T  is  defined.  The  macroscopic  quantities  are  given  by 
Euler  equations. 

3.1.1  General  non-equilibrium  situations 

The  most  general  non-equilibrium  situation  corresponds  to  the  equation  JRR  =  0 ,  that  is  the 

collisions  I  include  the  T  (elastic)  and  R  collisions  and  the  collisions  II  include  all  V  collisions. 
Thus,  we  have 

gi  exp (-sJkT) 


r  0 

fi  =nL 


f  m 


2nkT 


exp 


r  i\ 

mu 


v  2  kTj 


Qk 


(10) 


The  corresponding  macroscopic  equations  are  the  Euler  equations  and  the  relaxation  equations 
giving  the  vibrational  populations  nj  ,  that  is 


dn.  d.n.  V 


dt 


—  +  - 


— — =  V  —n:n:a 

&  ir,j,k,i 


k,l 


(ID 
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where  V  is  the  macroscopic  velocity,  a.J  and  a‘{d  the  collision  rate  coefficients  corresponding  to 
the  transitions  i,  j  I  k,l ,  independent  of  the  populations. 

3.1.2  Particular  non-equilibrium  situations 

Now,  if  particular  highly  probable  collisions  with  vibrational  exchanges  are  included  in  collisions 
I,  the  vibrational  populations  take  specific  forms.  Thus,  three  important  examples  are  briefly 
presented  hereafter. 

•  TV  collisions  are  dominant  and  are  included  in  collisions  I;  we  have: 

Jjy  ^ 

In  this  case,  the  solution  is  given  by  the  relation  (9),  which  means  that  the  TV  collisions 
are  sufficient  to  establish  an  equilibrium  regime. 

•  VV  collisions  are  dominant  [15]  and  are  included  in  collisions  I;  in  this  case,  we  find  a 
non-equilibrium  distribution  (10),  with 

exp  ( -  c:  jkT  +  Kiv  \ 

n.  =n— — y—  '  - (Treanor  distribution)  (12) 

‘‘  XexP  (-£jkT  +  Kiv) 

iv 

The  macroscopic  parameter  K  is  given  by  the  following  unique  relaxation  equation: 

^  =  [exp  (-/f )  -1  ]  X  <+i  Vi  ’ 

where  Ivis  a  mean  quantum  number  defined  as  nlv  -  ^  n.  /  ,  and  where,  here,  only 

iv 

monoquantum  transitions  are  assumed  to  take  place. 

This  distribution  presents  a  minimum  for  a  particular  quantum  number,  so  that 
“population  inversions”  are  possible  ,  especially  in  gas  mixtures. 


Resonant  collisions  (iv  +  jv  □  jv  +/v)  are  dominant  and  are  included  in  collisions  I.  In 
this  case,  we  find  a  Boltzmann  distribution  at  a  “vibrational  temperature”  Tv  ,  that  is 

exP  (-£JkTv) 


ni  =n 


ZexP  (-£JkTv) 


(13) 


This  temperature  is  also  given  by  a  unique  relaxation  equation  giving  the  mean  vibrational  energy  Ev . 

In  order  to  complete  the  previous  equations,  a  physical  oscillator  model  is  required.  For  example,  the 
choice  of  the  simple  harmonic  oscillator  model,  leading  to  the  Landau-Teller  relaxation  equation 
simplifies  the  problems  but  may  mask  detailed  (and  essential)  features  of  the  relaxation. . . 

On  the  other  hand,  before  choosing  a  global  model  (Boltzmann,  Treanor,. . .),  it  is  essential  to  determine  its 
domain  of  validity:  thus,  for  example  at  low  temperature,  the  VV  collisions  are  more  probable  for  the  low 
levels  and  the  TV  collisions  more  probable  for  the  highest  levels  at  high  temperature:  An  example  [16]  is 
presented  in  Fig. 3. 
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<;10(W),  a£~l  (TV )  (cmV) 


0  10  20  30  /v 


Figure  3:  TV  and  VV  collision  rates  (N2  -N2  ,  T=1000K) 


From  a  general  point  of  view,  the  macroscopic  parameters  appearing  in  ft0  are  given  by  the 

conservation  equations  including  Euler  and  relaxation  equations  corresponding  to  the  Fredholm 
alternative  of  the  first  order  equation  of  the  system  (8),  that  is 

\^r-OC,dv=0  , 

V 

where  C,  are  the  collisional  invariants  of  collisions  I. 

3.2  First  order  solutions.  Transport  Terms 


3.2.1  WNE  Case 


This  is  a  relatively  classical  case.  We  start  from  the  equilibrium  zero  order  solution  (9)  and  we 

7  r  0  t  r  0 

have  to  solve  the  equation  — —  =  J'TKV  .  Taking  into  account  the  form  of  the  known  terms  — — 

dt  dt 

and  various  constraints  related  to  the  definition  of  macroscopic  parameters  with  ft° ,  we  can  write 
the  perturbation  (pj  in  the  following  form: 


„  1  8T  n  dV  0  „  d.V 

C 'P- ;  =  A - .u  +  B;  - :  uu+  D;  - 

T  or  dr  or 


(14) 


where  Xt  =  A,  Bn  D  are  unknown  scalar  functions  of  r,  t,  u,  £i  ,  st  .  Classically,  each  term  Xt 
is  expanded  in  Sonine-Wang-Chang-Uhlenbeck  polynomials  [17],  that  is 


o mpnpp 
^mnp^T  1  R  rV 


(15) 


m,n,p 


with  x  =  a,  b,  d .  In  the  expansion  (15),  only  zero  and  first  order  terms  are  retained,  i.  e.  x000  , 
a‘| oo ,  v0|{| ,  x001 ;  each  index  corresponds  to  the  order  of  expansion  respectively  for  translation, 
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rotation  and  vibration.  These  terms  are  determined  from  particular  equations  involving  the 

df° 

corresponding  terms  of  — — .  Because  of  the  properties  of  the  orthogonal  expansion  (15),  they 

dt 

can  be  expressed  as  functions  of  “collisional  integrals”,  defined  as  follows: 


(...) 


kT 


X  ^  J  exp (-r2)r3 (...)!■  j  do.  dr  ,  (16) 


\7lTYl  J  i,j,k,l  ft  Q  y 


where  the  populations  n. ,  nj  are  in  equilibrium  and  where  I-'j  is  the  differential  cross-section  of 


the  transition  i,  j  —> k,l,  y  = 


f  m  \/2 


2  kT 


g  a  dimensionless  relative  velocity  and  Q  the  diffusion 


solid  angle.  These  integrals  involve  all  types  of  collision. 

At  the  macroscopic  level,  the  Navier-Stokes  equations  govern  the  evolution  of  macroscopic 
quantities  n,  V,  T  and  the  corresponding  transport  terms  are  the  following: 


8V  d.V 

Stress  tensor:  P  =  pI  +  P  ,  with  :  P  =  -2 u - n - 1 

dr  dr 


Heat  flux: 


q=qT+qR+qv  =~* 


dT_ 

dr 


(17) 


where  the  viscosity  and  conductivity  coefficients  //,  77,  X  are  also  known  functions  of  collisional 
integrals  of  type  (16).  More  details  and  approximate  expressions  may  be  found  in  Ref.24  .  Thus, 
we  can  see  that  pUrTp  and  //  U  rv  p .  It  is  important  to  note  that  a  “weak”  rotational  and 
vibrational  non-equilibrium  appears  at  first  order  of  the  expansion  of  the  distribution  function: 
since  generally  rv  □  rR,  we  find  for  the  vibrational  non-equilibrium: 

dV 


Ey  □  Ty 


dr 


(18) 


Temperatures,  specific  to  each  energy  mode  Tt,  Tr,  Tv  may  thus  be  defined  at  first  order  but  this 
is  independent  of  the  definition  of  a  single  “average”  temperature  T,  common  to  the  three  modes 
and  given  by  Navier-Stokes  equations.  It  may  be  also  noted  that  the  difference  between  Tt  and 

Tr  is  of  the  order  of  tr  ,  so  that  it  is  generally  possible  to  define  a  common  temperature  Tm , 


which  is  of  course  not  the  case  for  Tt  and  Tv . 


3.2.2  SNE  Case 

We  start  from  the  zero  order  non-equilibrium  solution  (10),  with  vibrational  population  given  by 
Eq.  (13)  in  which  the  dominant  collisions  include  the  vibrational  resonant  transitions:  this  enables 
us  to  preserve  the  notion  of  vibrational  temperature  Tv .  Thus,  the  first  order  equation  is  the 
following: 

^  =  jlTR+j0y  (19) 
dt 

df° 

Taking  into  account  the  zero  order  terms  — — Jy,  we  can  write  the  perturbation  <pi  in  the 

dt 

following  form: 
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.  1  dT  D  8V  0  ^  d.V  1  dTy 

a),  =  A - .u  +  B- - :  uu+  D, - 1-  F - -.u  +  G, 

T  or  or  or  Tw  or 


(20) 


where,  as  above,  Xt  =  An  Bn  Df,  F,  G;  are  unknown  scalar  functions  of  r,  t,  u,  et ,  et  . 

The  last  two  terms  in  Eq.(20)  arise  from  the  vibrational  non-equilibrium  already  present  at  zero 
order. 

Expanding  Xt  as  above  (15),  we  can  also  express  the  coefficients  x  as  functions  of  collisional 
integrals:  however,  as  expected  from  Eq.(19),  these  coefficients  (except  g)  depend  on  integrals 
including  only  the  collisions  I,  that  is  here  the  TR  collisions,  noted  (...)TR.  The  coefficients  g 

directly  depend  on  the  vibrational  relaxation,  because  of  the  7°  term. 

The  macroscopic  conservation  equations  are  also  the  Navier-Stokes  equations,  completed  by 
relaxation  equation  giving  ev  or  Tv  ,  that  is 


dE, 


n- 


v  !  d4v 

dt  dr 


=2>JPv0+4)a- 


Of  course,  a  physical  oscillator  model  is  necessary. 

Moreover,  the  transport  terms  in  the  Navier-Stokes  equations  become 


(21) 


Stress  tensor:  P  =  pl-2ju 


dV 

dr 


dV 
V - 

v  dr 

Heat  fluxes  :  q  =  qT  +qR  +qv  ,  with 
_  dT  dTv 

qT  Ay  Ayy 

dr  dr 

_  dT  dTv 

Qr  ^ R  ^  A RV  ^ 

dr  dr 

_  dT  dTv 

qv  A'tr  ^  /V  o 

dr  dr 


Pr 


(22) 


(23) 


Thus,  in  the  stress  tensor  (22)  we  find  the  “usual”  terms  and  a  “relaxation  pressure”  term  pR 

(related  to  7°).  In  the  heat  flux  expressions  (23),  we  find  cross  terms  depending  on  both 

temperature  gradients.  However,  as  already  discussed,  the  transport  coefficients 
p,  77,  A... include  only  TR  collisions  (and  resonant  V  collisions).  The  bulk  viscosity  rj  for 

example  depends  only  on  zR . 

Thus,  starting  from  zero  order  solutions,  either  in  equilibrium  or  not,  we  find  structural 
differences  in  the  description  of  non-equilibrium  situations  at  first  order.  This  is  summarized  in 
Table  I. 
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Method 

WNE 

SNE 

Zero  order 

J°  =0 

J  TRV  V7 

TRV  equilibrium,  with  single 
temperature  T. 

Euler  equations 

j°  =o 

JTr  V 

TR  equilibrium,  with 
common  temperature  T. 
Vibrational  non-equilibrium,  with 
temperature  Tv. 

Euler  equations,  with 
vib.  relaxation  equation 

First  order 

df-  _  ,i 

J  TRV 

dt 

Weak  RV  non-equilibrium. 

Navier-Stokes  equations 

df‘°  =  J1  +  jo 
j  TR-r  Jy 

dt 

Weak  R  non-equilibrium. 
Navier-Stokes  equations,  with 
vib.relaxation  equation 

Transport 

Transport  terms  //,  i),  A 
depending  on  TRY  collisions. 
Heat  fluxes  depending  on  the 
gradient  of  T 

Transport  terms  //,  77,  Aj ,  AR ,  Ay 

depending  on  TR  collisions 
Heat  fluxes  depending  on  the 
gradients  of  T  and  Tv 

Table  I 

Comparison  of  WNE  and  SNE  methods  for  vibrationally  relaxing  gases 


It  is  thus  obvious  that  there  is  no  matching  between  both  first  order  solutions  when  Tv  — »  T .  That 
is  why,  a  generalized  Chapman-Enskog  method  (GCE)  is  proposed  hereafter. 

3.2.3  GCE  Method 


The  matching  between  WNE  and  SNE  cases  is  simply  realized  by  adding  a  Jxn  term  in  the 
linearized  equation  of  the  SNE  system  (8),  so  that  this  system  becomes: 

J°=0 


dt 


=  jl  +  j°  +  j1 

J  ,  T  J  u  T  J  jj 


(24) 


At  zero  order,  we  have  the  SNE  solution  covering,  as  previously  discussed,  the  non-equilibrium 
and  equilibrium  situations  (Euler  system  +  kinetic  and/or  relaxation  equations).  At  first  order,  far 
from  equilibrium  the  J),  term  is  negligible  and  we  find  the  SNE  solution,  but  close  to 


equilibrium  J°n  — >  0  and  the  first  order  solution  is  the  WNE  solution. 

The  first  order  equation  of  the  system  (24)  is  however  more  difficult  to  solve  than  in  the  WNE  or 
SNE  cases  because  of  the  structure  of  the  operator  J), .  Thus,  this  operator,  contrary  to  j)  is  not  a 
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self-adjoint  operator,  but  it  may  be  decomposed  into  two  operators;  among  them  the  dominant 
term  is  self-adjoint  and  the  other  may  be  neglected.  Finally,  in  the  present  case  of  vibrational 
non-equilibrium,  we  can  write  the  system  in  the  following  way: 


/  =0 
JTR  J 


df[_ 

dt 


=  jl  +J°+  JA1 

°  TR  '  u  V  '  u  V 


(25) 


As  previously  discussed,  the  zero  order  solution  of  the  system  (25)  is  the  SNE  solution  (13)  with 
Euler  and  relaxation  equations.  The  perturbation  (pt  is  given  by  Eq.(20)  and  we  have  the  same 

equations  for  the  coefficients  X j ,  but,  in  these  equations,  the  operator  J'IK  is  replaced  by 
JjR  +  Jyl .  Thus,  solving  the  equation  systems  giving  the  coefficients  x,  we  have  to  add  collisional 
integrals  to  collisional  integrals  (...)TR  similar  to  integrals  (16),  that  is: 


(26) 


(  kT  \/2  __ 

—  I 

V  f  exp(  y  )r  {..)L  ’j(TR)dQ.dy 

Knm)  iJtktl 

n  n 

_  Q,/ 

<->v 


\nm  J 


r 

f 

f  tV 

\  T 

i 

l  +  exp 

A  6* 

1 - 

X 

Z- 

V 

V 

v 

)  > 

n 

Z 

|  exp| 

l-r) 

ir3(-K 

cj  (V)dQdy 

Q.,/ 

,(27) 


Sk  +£,  ~8i  -£■ 

where  A  sv  —  — - - - - - - 11  ^  0  and  where  n,  ,  tij  are  vibrational  non-equilibrium 

kT 

populations. 

Thus,  in  integrals  the  non-equilbrium  appears  explicitly.  Navier-Stokes  equations,  of 

course,  remain  valid  as  well  as  relaxation  equations  similar  to  Eq.(21),  but  the  transport 
coefficients  such  as  /u,  t],  /(...depend  on  the  relaxation  explicitly.  For  example,  the  viscosity 
coefficient  /u  may  be  written  as  follows: 

—  Rtr  z  tr  ,  (28) 

1  +  K 

where  R  represents  the  ratio  of  two  integrals  similar  to  (26)  and  (27),  that  is  R  =  . 

(0)tr 

For  the  heat  transfer  terms,  we  find  expressions  similar  to  those  of  the  SNE  system  (23),  but, 
when  Tv  — ,  we  find  the  expressions  of  the  WNE  case  (17)  for  the  conductivity  coefficients. 

Approximate  expressions  may  be  obtained  and  the  GCE  method  can  be  extended  to  gas  mixtures: 
details  may  be  found  in  Ref.  10  .  An  example  of  vibrational  conductivity  is  represented  in  Fig.4. 
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0  1  2  3  T  (103£) 

Figure  4:  Vibrational  conductivity  of  nitrogen  in  a  mixture  (air) 

( A:TV  =3000K,B:TV  =1000K,C:TV  =800 K,D:TV  =500 K,E:TV  =T ) 


4.0  DISSOCIATING  GASES 


4.1  (WNE)v+(SNE)c  Case 

We  consider  dissociating  polyatomic  gases;  as  generally  the  characteristic  chemical  times  are 
larger  than  vibrational  relaxation  times  (Figs.l  and  2),  we  consider  a  WNE  case  for  vibration  and 
a  SNE  case  for  chemistry,  that  is  a  (WNE)v+(SNE)c  case.  For  a  molecular  component  p,  we  have 
therefore  the  following  system: 


J  TRVp  0 


dt 


=  JLu„  +  J, 


(29) 


'  TRVp 


Cp 


At  zero  order,  we  have  a  Maxwell-Boltzmann  distribution  (9)  and  corresponding  Euler  equations 
closed  by  kinetic  equations  of  the  following  type: 

dPp  .  d-PjY  .  Y*  f  TO  7 

=  wP=Lmp\JcPdvp  (3°) 


dt 


-  +  - 


dr 


Considering  first  only  a  dissociating  pure  gas,  we  have  two  species,  molecules  p  and  atoms  q. 
Neglecting  the  recombination  process,  we  simply  have 


w„ 

— —  =  -knn 


m 


Dtlp 


(31) 


where  kD  =kD(T)  is  the  Arrhenius  dissociation  rate  coefficient. 

At  first  order,  we  can  write  the  perturbation  <pjp  in  the  following  way: 


1  dT  dV  0  dV 

(p.  =  A. - .u  +B  - :u  u  +D - 1 -G  +L  u  <d 

Yip  ip  T  dr  p  ip  dr  p  p  ip  dr  ip  "  ” 


ip  p 


(32) 


Comparing  with  the  previous  expressions  for  cp.  ,  we  see  in  Eq.(32)  a  diffusion  term  L.  . 
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The  transport  terms  of  the  Navier-Stokes  equations  may  be  calculated  as  above  but  also  the 
dissociation  rate  constant  at  first  order;  thus  we  find 


kD  kD 


1- 


£ooi  + 


001 


d_V_ 

dr 


V Ey-En ^ 


kT 


(33) 


Similarly,  we  find  for  the  vibrational  energy  at  first  order 

~  (  ay  3 

Ey  =  Ey  1  +  goOl  ^001  “T 
V  dr 


(34) 


so  that  combining  (33)  and  (34),  we  obtain  an  expression  for  the  rate  constant  depending  on  the 
vibrational  non-equilibrium,  i.e. 


kD  kD 


1  + 


75  Vr 


-E  ^ 


kT 


(35) 


where  EVD  represents  the  vibrational  energy  loss  due  to  dissociation  (at  zero  order). 

We  need  a  physical  model  of  course  in  order  to  have  a  complete  expression:  details  may  be  found 
in  Ref.ll. 

4.2  (SNE)v+(SNE)c  Case 

We  consider  simultaneous  vibrational  and  chemical  non-equilibriums.  Thus,  we  have 


j°  =0 

JTRp  u 

df° 

_Ev_  —  j1  +J°  +J° 

J  TRp  '  J  Vp  '  J  Cp 


(36) 


At  zero  order,  we  have  simultaneously  a  vibrational  relaxation  and  a  chemical  production. 

This  case  has  been  widely  analysed  in  the  past  at  the  macroscopic  level  [19],  [20],  [21]: 
dissociation  rate  constants  are  modified  and,  depending  on  various  assumptions  and  models, 
expressions  for  kD  =  kD  (T,  Tv  )  have  been  obtained.  Similarly,  the  usual  vibrational  relaxation 
equation  is  modified  because  of  the  dissociation. 

However,  at  first  order,  no  important  interaction  takes  place,  since  the  first  order  terms  depend 
only  on  rotational  relaxation,  generally  negligible.  It  is  the  same  for  the  transport  terms  which  do 
not  depend  on  vibration  (see  previous  SNE  cases).  Further  discussions  may  be  found  in  Ref.24. 

kD(T,Tv) 

An  example  of  comparison  for  the  ratio  V  =  —  computed  from  different  models  is 

represented  in  Fig.5. 

4.3  Tentative  general  solution 

In  an  attempt  to  unify  the  calculation  methods  of  relaxing  and  reacting  gas  flows  at  the 
macroscopic  level,  we  suggest  to  use  the  Navier-Stokes  equations  closed  with  species 
conservation  equations  (kinetic  equations)  and  vibrational  relaxation  equations  taking  the 
chemical  reactions  into  account:  this  is  an  SNE  formulation,  but  the  reaction  rate  constants  are 
computed  with  the  WNE  model  [8], 

For  example,  in  the  simple  case  of  dissociating  pure  gases,  the  formula  for  kD  is  given  by  Eq.(35) 
and  the  relaxation  equation  may  be  written  as  follows,  if  we  use  the  harmonic  oscillator  model: 


dEv 

dt 


Ey  ~Ey 


+  (  Ey  Evd  j  kDnp 


(37) 
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The  term  EVD,  is  calculated  at  first  order,  like  .  We  can  note  that  realistic  values  for  EVD  lie  in 
the  range  0.3-0.5  ED  [13]. 

In  this  framework,  the  transport  terms  may  be  used  in  their  GCE  formulation  which  takes  into 
account  the  vibrational  relaxation.  From  a  general  point  of  view,  it  is  also  recommended  to  use 
these  terms  in  their  dimensionless  form,  less  sensitive  to  non-equilibrium  effects. 

Extensions  to  dissociating-recombining  case  and  to  gas  mixtures  including  atomic  and  molecular 
components  as  well  as  numerous  reactions  have  been  carried  out.  A  few  examples  are  given 
below. 


LogV(T,Tv) 


Sc 

T=2.\tfK 

0  4  8  12  16  20  Tv  (lO3^) 

Figure  5:  Vibration-dissociation  factor  V  for  nitrogen,  T=2.104  K 

A:  (SNE)C+(WNE)V  Case  :  Non-preferential  anharmonic  oscillator  model 

B:  (SNE)C+(WNE)C  Case  :  Preferential  model,  -U=ED/6k  [20] 

C  :  Semi-empirical  model  [21] 


5.0  APPLICATIONS  TO  HYPERSONIC  FLOWS 
5.1  Shock  Waves 

We  consider  the  flow  behind  a  straight  strong  shock  wave  moving  in  quiescent  air,  with  the  following 
conditions:  Ms  =25,  p0  =8.5 Pa,  =205 K .  A  classical  chemical  model  including  17  reactions  is 

adopted  and  the  classical  Landau-Teller  equations  for  the  molecular  components  are  used;  however, 
formulas  of  the  type  (35)  are  used  for  the  chemical  rate  constants  involving  dissociation  and  exchange 
reactions  and  the  relaxation  equations  include  energy  loss  (gain)  due  to  chemical  reactions  (cf  Eq.37). 
One-dimensional  Euler  equations  are  also  used  [12], 

An  example  of  variation  of  a  dissociation  rate  constant  ( N2 )  behind  the  shock  is  shown  in  Fig.6  and 
compared  to  the  corresponding  Arrhenius  rate  constant.  Thus,  strong  differences  may  be  pointed  out  due 

to  the  vibrational  relaxation  in  the  region  where  Tv  (  N2 )  <  T . 
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The  evolution  of  the  vibrational  temperatures  is  represented  in  Fig.7  where  we  can  see  a 
minimum  for  Tv  (02)  which  does  not  exist  without  vibration-chemistry  interaction. 


In  order  to  give  an  idea  about  the  order  of  magnitude  of  the  “relaxation  pressure”  pR,  the 
evolution  of  is  represe 
the  order  of  a  few  per  cent. 


evolution  of  is  represented  in  Fig.  8:  thus,  close  to  the  shock  wave,  we  see  that  this  ratio  is  of 


kD  (m’ .mole  'js  1 ) 


0  0,4  0,8  1,2  1,6  2,0  x  (cm) 


Figure  6:  Dissociation  rate  constants  of  nitrogen  behind  a  shock  wave  in  air 

Ms  =25,  p0=8,5Pa,  T0=205K 


A:  kD  Arrhenius,  B:  kD  with  interaction,  C:  kD  empirical 


T  (104AT) 

25 

20 

15 

10 

5 

0 

0  2  4  6  8  10  x  (cm) 

Figure  7:  Spatial  variation  of  temperatures  behind  a  shock  wave 
in  air  (Conditions  of  Fig.6) 
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Pr/P 


Figure  8:  Relaxation  pressure  behind  a  shock  wave  in  air 

(Conditions  of  Fig.6) 


5.2  Boundary  Layers 

We  consider  the  boundary  layer  in  vibrational  and  chemical  non-equilibrium  developing  at  the 
end-wall  of  a  shock  tube  after  the  reflection  of  the  incident  shock  wave.  The  reflected  shock  wave 
propagates  in  the  non-equilibrium  medium  generated  by  the  incident  shock  wave:  an  example  is 
represented  inFig.9,  [13]. 

What  is  of  interest  here  is  the  behaviour  of  the  rate  constants  in  the  end-wall  boundary  layer;  thus, 
the  time  evolution  of  the  vibrational  non-equilibrium  across  the  boundary  layer  is  represented  in 
Fig.  10  in  which  we  can  see  a  vibrational  freezing  zone  ( Tv  >T)  appearing  some  time  after  the 

reflection,  due  to  the  migration  towards  the  wall  of  high  energy  molecules.  The  result  is  an 
increase  of  the  dissociation  rate  constant  in  this  zone  in  spite  of  the  decrease  of  T  (Fig.  11).  The 
phenomenon  is  accentuated  when  the  wall  is  non-catalytic. 


0  10  20  30  40  50  x  (cm) 

Figure  9:  Temperature  distributions  behind  a  reflected  shock  wave 

(Nitrogen,  Ms=12,  p0=200Pa) 

Time  after  reflection  :1  :15ps,  2  :25ps,  3  :40ps,  4  :400ps,  5  :500ps,  6  :600ps,  7  :700ps 
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0  1  2  3  4  5 


V 


Figure  10:  Vibrational  non-equilibrium  distribution  across  the  boundary  layer 

(Conditions  of  Fig9,  y/  :  Von  Mises  transverse  coordinate), 

Time  after  reflection  :  A  :25ps,  B  :500ps,  C  :  900ps 


Full  line  :  Catalytic  wall.  Dotted  line  :  Non-catalytic  wall 


Log(kD/kD) 

1,5 

1,0 

0,5 

0,0 

-0,5 

-1,0 

0  1  2  3  4  5  y/ 


Figure  11:  Dissociation  rate  constants  across  the  boundary  layer 

(Conditions  and  notation  of  Fig. 9) 


5.3  Flows  around  Bodies 

Comparisons  of  results  obtained  from  Navier-Stokes  computations  using  the  present  model  discussed 
above  with  experimental  results  carried  out  in  shock  tunnel  and  ballistic  range  are  presented  in  Figs.  12  and 
13.  In  fact,  they  refer  to  the  shock  stand-off  distance  over  hemisphere-cylinder  bodies:  this  parameter  is 
sensitive  to  the  non-equilibrium  effects,  as  it  may  be  seen  in  the  figures  [14].  The  good  agreement 
observed  between  experimental  and  simulated  values  gives  support  to  the  reliability  of  the  models  used. 
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Q  (Reaction  rate  parameter) 

Figure  12:  Measured  (•)  [22]  and  computed 
stand-off  distances  in  carbon  dioxide  flow  (shock  tunnel) 


2000  2500  3000  3500  4000 


Flight  velocity,  m/s 

Figure  13:  Measured  [23]  and  computed  shock  stand-off  distances 

pR  =  2.  NT4  kg  /  m2  ,  A  :  Experiments  (gun  tunnel) 
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6.0  CONCLUSIONS 

Much  work  remains  to  be  done  in  the  field  of  kinetic  theory  of  reactive  molecular  gases:  more  elaborate 
models,  STS  models,  inelastic  and  reactive  cross-sections,  electronic  transitions,  direct  simulations, 
complete  GCE  method,  ....  The  results  presented  above,  however,  prove  a  relative  maturation  of  ideas  and 
methods,  so  that  reliable  computations  of  non-equilibrium  flows  may  be  carried  out. 
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